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We develop a theory to compute and interpret the photonic band structure of a periodic array of metalhc 
hehces for the first time. Interesting features of band structure include the ingenuous longitudinal and 
circularly polarized eigenmodes, the wide polarization gap [Science 325, 1513 (2009)], and the helical 
symmetry guarantees the existence of negative group velocity bands at both sides of the polarization 
gap and band crossings pinned at the zone boundary with fixed frequencies. A direct proof of negative 
refraction via a chiral route [Science 306, 1353 (2004)] is achieved for the first time by measuring Goos- 
hanchen shift through a slab of three dimensional bona fide helix metamaterial. 

PACS numbers: 78.67.Pt, 42.70.Qs, 42.25.Ja 



With a pitch angle accounting for mirror asymmetry, a 
helix is a prototype element of chiral media. As a pi- 
oneer study, the electromagnetic(E]V[) activity from ran- 
domly dispersed metallic helices was observed in 1914 
fit]. After the prediction on negative refraction via chiral 
route [2], chiral metamaterials have been arousing more 
and more attentions recently on negative refraction index 
IS Mol, s trong optical activity [7-9] and circular dichroism 
I1Q|-[13Ii as well. In the most of previous studies, chiral 
metamaterials are typically realized as one or two layers of 
discrete chiral resonators and the negative refractive index 
is deduced indirectly from retrieval procedures as- 
suming the applicability of effective medium theory 11141 - 
[iS]. Recently, the Gold helix metamaterial with a few 
pitches along the helical axis are successfully fabricated in 
THz regime with a wide polarization gap measured exper- 
imentally [19], which also raised a question whether the 
effective medium theory is enough for the description of 
chiral metamaterials. To the best of our knowledge, most 
previous studies adopt effective medium description and 
well-defined EM parameters for theoretical analysis, and 
empirically interpret the exotic properties in a numerical 
stage. While there is no any "first-principle" investigation 
starting from the real structure of chiral metamaterial. Also 
it is worth noting that in all previous studies, negative re- 
fraction indices, retrieved from the transmission and reflec- 
tion spectra under normal incidence, do not give any direct 
proof for negative refraction. The negative refraction via 
chiral route hasn't really achieved before. 

In this paper, we theoretically and experimentally inves- 
tigate the properties of a three-dimensional metamaterial 
made with a square array of metallic helices. We develop a 
photonic band theory for helix metamaterial by combining 
IVLultiple Scattering Theory (IVLST) [20, 21] with Sensiper's 
solution for a single helix ll22|] . The theory enables us 
to identify the longitudinal and circularly polarized eigen- 
modes, the formation of the wide polarization gap, negative 
group velocity bands at both sides of the gap and related 
fine features. We also demonstrate a proof of principle 



experiment by measuring negative Goos-hanchen shift to 
directly verify the negative refraction from helix metama- 
terial for the first time. 

Figure [T] shows a photograph of a sample of the three- 
dimensional metamaterial comprising a square array of 
right-handed (RH) helices. Oriented along the z axis, the 
metallic helices form a square array in the xy plane with 
a lattice constant of d = 11 mm. Figure [TJa) shows the 
schematic picture of a single helix unit, which has a pitch 
of p — 4.4 mm, radius a = 3.3 mm, and diameter of 
metallic wires 5 = 0.8 mm. We define a pitch angle by 
cot ^ = 27va/p. We note that a helix comes back to itself 
after being translated by a distance of Az and being rotated 
simultaneously by a angle of 27rAz/p (for RH helix) or 
—27rAz/p [for left-handed (LH) helix], and thus physical 
entities associated with the helix should satisfy the helical 
symmetry condition 
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with +/- sign for the RH/LH helix respectively SIH]. As 
the helical system is also periodic along the helical axis, the 
field components for an RH helix system can be expanded 
by functions of the form 
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FIG. 1 : (a) The schematic picture of a metallic helix unit and (b) 
photo of a square array of metallic helices. 




FIG. 2: The photonic band structures of the hehx crystal along the helix axis for three lattice constants (a) d = 20 mm, (b) = 11 mm, 
(c) d = 8 nmi. The pitch p = 4.4 mm, the radius a = 3.3 nrni, and diameter of metallic wires ^ = 0.8 mm. 



where is the Bloch wavevector along z axis. The angu- 
lar term should be e*^"^ instead if the helix is LH. The radial 
function Fn (p) obeys the Helmholtz differential equation, 
and can be expressed in terms of modified Bessel func- 
tions In and We follow the Sensiper approach [22] to 
impose the assumption about uniformly distributed surface 
current flow along the metal wires. Under that assumption, 
the boundary continuity conditions require that the local 
electric field on metal wires must be perpendicular to the 
line of metal wire. Thus we derive an eigenvalue equation 
of helix array by MST as 
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Si{r) = Ki{rRq)e'^'f'^e'^'-^ is a lattice sum 

running over the nodes {Rq^(j)q) of the square lattice in 
cylindrical coordinates with ki being the transverse com- 
ponent of the wavevector k in the air, and /^(x), K'^[x) 
satisfy to /^(x) = dln{x)/dx, K'^(x) = dKn{x)/dx. 

The photonic band structures, computed with the com- 
bined technique above stated, afford us an intuitive under- 
standing about optics of helix metamaterials. In the helical 
structure, there is always a7r/2or— 7r/2 phase difference 
between the radial and the angular components for both the 
electric field and magnetic field, implying that the eigen- 
modes are ingenuous left-handed or right handed circularly 
polarized (LCP or RCP). This can be checked by exam- 
ining the field solutions written in the Sensiper form[22] 
, or by examining the solution from commercial numeri- 
cal solver CST. Figures Oa), Ob) and^c) show the band 
structure of the helix metamaterials along the helix axis for 



three lattice constants d = 20 mm, 1 1 mm and 8 mm re- 
spectively. We label the eigenmodes by its dominant term 
in Eq. (O. For example, the (-1,S) modes (blue stars in 
Fig. [2]) have n = — 1 term as the dominant term and "S" 
stands for a "slow mode" below the light line, and we use 
the subscript "F" for a mode inside the light cone. The low- 
est frequency branch (red circles), labeled as (0) in Fig. [2l 
has a strong longitudinal component and this mode picks 
up a LH character as k^ increases. It goes to a longitudinal 
mode with a finite frequency at Brillouin zone (BZ) center 
(kz = 0). Both the electric field and magnetic field are 
essentially parallel to the helical axis. It is evident from 
Fig. [21 that the inter-helix coupling pushes the longitudi- 
nal mode to higher frequencies. In the limit d oo, it 
goes to zero frequency. The polarization of an eigenmode 
is analyzed by the ratios \{E,)\ / \{E,)l \{E,)\ / \{Ey)\ 
and AR = (E^) k^/ {iEy) | A;^ |, where the spatial average, 
(•••), is taken inside a unit cell. Figure [3] clearly indicate 
that the polarization of n = 0, ±1 eigenmodes are longi- 
tudinal and/or LCP, RCP exclusively. Eigenmode analysis 
and numerical transmission simulation employing a finite 
thickness slab showed that the n = — 1 modes with posi- 
tive group velocity (blue stars in Fig. O couple to incident 
plane wave with opposite handedness as the helix, while 
the n = +1 modes with positive group velocity (green 
squares in Fig. O couple to incident waves of the same 
handedness as the helix. In general, a mode couples to in- 
cident wave of the same (opposite) handedness as the helix 
if n-k,> 0(< 0). 

An important feature of the band structure is the exis- 
tence of a wide polarization gap (shaded in grey in Fig. O 
that only allows incident waves of opposite handedness to 
pass through. As such, a RH helix array has a right-handed 
circularly polarized (RCP) gap. The gap grows wider for a 
higher helix filling ratio (smaller d). Such wide polariza- 
tion gaps have been experimentally demonstrated for a thin 
slab of gold helices in IR frequencies Ji^ll^]. The lower 
edge of the polarization gap is pinned at the frequency at 
which the (+1,S) mode attains zero group velocity. The 
(+1,S) mode is the result of the hybridization between the 



3 



5x10 ■ 
4x10^ 
3x10' 
2x10' 
1x10' 


-0, 

3 
2 
1 

0- 
-1- 



- |<Ez>|/|<Ex>| 
-|<Ez>|/|<Ey>| 



(0) Modes 



(0) Modes 



0.00 



0.25 



(+1,S) Modes 



'-^ |<Ez>|/|<Ex>| J 
|<Ez>|/|<Ey>| T 

-■-Re(AR) / - 
Im(AR) / 



|<Ez>|/|<Ex>| 
- |<Ez>|/|<Ey>| 
0-|«»*********«*«**««>«*«*««*««V«**«***« 

-■- Re(AR) 
-#- Im(AR) 



-0.25 



0.00 



0.25 



-1 
-0.50 



(+1,F) Modes 



-0.25 



0.00 



0.25 



\ |<Ez>|/|<Ex>| 
V |<Ez>|/|<Ey>| 
V -■- Re(AR) 
%-*-lm(AR) 



(-1,S) Modes 



|<Ez>|/|<Ex>| 
■^|<Ez>|/|<Ey>| 

-■- Re(AR) 
-#- Im(AR) 



-0.25 0.00 
k a/cot 



-0.25 0.00 0.25 



FIG. 3: Eigenmode analysis by calculating the ratios AR = 

{Ea:)kz/{iEy)\kzl \{E,)\ / \{Ea:)\ and \{E,)\ / \{Ey)\. (a) and 
(b) for n = modes; (c) (+1,S) modes; (d) (+1,F) modes; (e) 
(-1,S) modes; (f) (-1,F) modes. 



free photons (riding on the light line) and a mode guided on 
the helix which is strongly back- scattered by the periodic- 
ity. As the (+1,S) mode is guided on the helices, the fre- 
quency is only moderately affected by the inter-helix cou- 
pling. The upper edge of the polarization gap is pinned 
at the frequency at which the (+1,F) mode has zero group 
velocity. The (+1,F) mode is not a guided mode on a sin- 
gle helix. It exists as an eigenmode only in the periodic 
array and is sensitive to the volume available between the 
helices and the mode is squeezed to higher frequencies at 
a higher helix filling ratio, leading to a much bigger polar- 
ization gap for smaller values of the lattice constant (d). In 
order words, the lower edge of the gap is primarily deter- 
mined by the helix parameters (a and p) along the helical 
axis, while the upper edge is primarily determined by the 
structural parameter perpendicular to the helical axis (lat- 
tice constant d). 

Another interesting feature of the band structure is the 
emergence of the negative group velocity bands at both 
sides of the polarization gap, which is different from the 
previous notions that the negative refraction only happens 
above the resonant gap. Both the high frequency (+1 ,F) and 
the lower frequency (+1,S) branch exhibit negative group 
velocities. The (+1,S) branch has negative group velocity 
after reaching a maximum frequency that pins the lower 
edge of the polarization gap. Concomitant with the nega- 
tive refraction bands in slow mode, one can see that there 
are no band gap but band crossings at the BZ boundary 
(k = 7v/p) and the degenerate modes are pinned at fixed 
frequencies that are nearly independent of lattice constant 
p. We note that according to Eq. [2]Bloch modes that differ 
by Ak = 27r/p will be orthogonal because of orthogo- 
nality of the angular phase factor. And the (0) branch de- 
generates with the (+1,S) branch at zone center, giving rise 
to a negative group velocity band below the polarization 



gap. This is essentially a consequence of backward waves 
guided on metallic helix satisfying to helical symmetry. As 
a proof, a comparative study shows that the negative dis- 
persion band disappears as well as other unique features, 
such as the longitudinal mode, band crossing etc. when the 
helix metamaterial is cut in a "discrete" form by inserting 
an 0.4mm air gap at each period(p). 

We performed negative refraction measurements inside 
an anechoic chamber through a slab of the helix metama- 
terial with the aforementioned geometric parameters [the 
band structure is shown in Fig. Ob)]. Helix metamate- 
rial samples are fabricated by periodically embedding the 
clockwise metallic helices in a polyurethane foam slab. 
The polyurethane foam slab is lossless with e ^ 1. The 
sample slab contains 15x11 metallic helices, each having 
140 periods along helical axis (z axis). 

Computed equi-frequency surface (EFS) for the (+1,S) 
branch [see Fig. H^b)] and (+1,F) branch (not shown) 
demonstrate that negative refraction can be achieved at 
both sides of the polarization gap. Here we try to real- 
ize the negative refraction at the lower edge of the gap, 
which is not found or predicted in other systems before. 
As the (±1,S) modes lie below the light line, we can excite 
the (±1,S) modes by prism coupling techniques to pene- 
trate a Gaussian beam into the sample slab [see Fig. HJa)], 
and estimate the refractive angle quantitatively by mea- 
suring the Goos-Hanchen shift of the beam. Figure Ufa) 
illustrates the schematic configuration of our experimen- 
tal setup. Two isosceles right-angled triangular alumina 
prisms (Sr = 8.9) are placed so that they touch the sample 
slab at both sides and a Gaussian beam is normally inci- 
dent in xz plane to the air-prism interface from a linearly 
polarized horn emitter (operating at 8.2-12.4 GHz with a 
gain factor of 24.6 dB), ensuring an incident angle of 45 ° 
from alumina to metamaterial. The local field intensity is 
measured by the LCP/RCP horn receiver as a function of 
the horizontal position in a precision of 1 mm per step. The 
Goos-Hanchen shift of the outgoing beam is quantified by 
measuring the peak position at the interface of prism. The 
coordinate origin in the horizontal position is aligned to 
the position of the horn emitter, marked by the dashed ver- 
tical line in Fig. (H^a). The negative refraction is measured 
from 9.18 GHz to 9.48 GHz with a refraction angle from 
— 17.44 ° to —50.11 ° , which is in good agreement with 
the computed equi-frequency surface (EFS) shown in Fig. 
lUb). The solid line in Fig. |4td) presents the spatial pro- 
file of local field intensity measured at 9.41 GHz. A peak 
value is measured at the horizontal position of —24 cm, 
corresponding to a refraction angle of —46.45 ° , roughly 
equal to —45.82 ° estimated by EFS analysis. Thus neg- 
ative refraction below the polarization gap is verified ex- 
perimentally. Negative refraction is also qualitatively ver- 
ified by finite-difference-in-time-domain (FDTD) simula- 
tions shown in Fig. Hfc). We note from Fig. Hfb) that in the 
frequency range of our interest, the dispersion is not only 
negative along k^, but it is also negative along k^ due to 
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FIG. 4: (a) Schematic of the experimental setup for negative refraction, (b) Computed equi-frequency surface (EFS)for the (+1,S) 
branch. The arrow refers to the direction of the refracted waves at 9.41 GHz under an incident angle at 45 ° . (c) FDTD simulations 
of negative refraction. In calculation, an RCP Gaussian beam with a frequency f = 9.41 GHz is normally impinging on the left prism. 
The chiral sample has 4x5 helices in the xy plane and 45 periods along the z axis. The direction of energy flow is indicated by the 
arrow, (d) Measured electric field intensity as a function of the horizontal position of the circularly polarized horn receiver, with (red) 
or without (black) the chiral sample and alumina prisms at 9.41 GHz. The two curves are normalized such that the magnitude of both 
peaks is unity. The blue dashed line at horizontal position of —11 cm refers to Goos-Hanchen shift with respect to ° refracted angle. 



Bragg scattering. 

In summary, we witness the mathematical beauty of 
wave propagation in metallic helix metamaterial by analyt- 
ical resolved photonic band structure. There are negative 
bands both above and below the polarization gap. Nega- 
tive refraction on the low frequency branch is demonstrated 
directly for the first time by a Goos-hanchen shift exper- 
iment. We note that the optics of helix metamaterial is 
governed by the collective selection on the guided modes 
on helices satisfying to helical symmetry and (evanescent) 
Bragg scatterings within helix lattice and in passing there 
is no easy way to describe the phenomena using effec- 
tive medium parameters. This work was supported by 
NSFC (No. 10974144, 60674778), HK ROC grant 600308, 
the National 863 Program of China (No.2006AA03Z407), 
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